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The Faddeev equation for three-body scattering at arbitrary energies is formulated in momen- 
tum space and directly solved in terms of momentum vectors without employing a partial wave 
decomposition. In its simplest form the Faddeev equation for identical bosons, which we are using, 
is a three-dimensional integral equation in five variables, magnitudes of relative momenta and an- 
gles. This equation is solved through Fade summation. Based on a Malfiiet-Tjon-type potential, 
the numerical feasibility and stability of the algorithm for solving the Faddeev equation is demon- 
strated. Special attention is given to the selection of independent variables and the treatment of 
three-body break-up singularities with a spline based method. The elastic differential cross section, 
semi-exclusive d(N,N') cross sections and total cross sections of both elastic and breakup processes 
in the intermediate energy range up to about 1 GeV are calculated and the convergence of the 
multiple scattering series is investigated in every case. In general a truncation in the first or second 
order in the two-body t-matrix is quite insufficient. 

PACS numbers: 21.45-l-v 



I. INTRODUCTION 



During the last two decades calculations of nucleon-deuteron scattering experienced large improvements and refine- 
ments. Here different techniques have been applied, Faddeev calculations in configuration space Q and momentum 
space , and variational calculations based on a hypcrspherical harmonic expansion |i (J| . It is fair to say that below 
about 200 MeV projectile energy the momentum space Faddeev equations for three-nucleon scattering can now be 
solved with high accuracy for the most modern two- and three-nuclcon forces. A summary for these achievements can 
be found in Ref. 0, S S S 13 • The approach described there is based on using angular momentum cigenstatcs for 
the two- and three-body systems. This partial wave decomposition replaces the continuous angle variables by discrete 
orbital angular momentum quantum numbers, and thus reduces the number of continuous variables to be discretized 
in a numerical treatment. For low projectile energies the procedure of considering orbital angular momentum com- 
ponents appears physically justified due to arguments related to the centrifugal barrier and the short range of the 
nuclear force. However, the algebraic and algorithmic steps to be carried out in a partial wave decomposition can be 
quite involved when solving the Faddeev equations. If one considers three-nucleon scattering at a few hundred MeV 
projectile energy, the number of partial waves needed to achieve convergence proliferates, and limitations with respect 
to computational feasibility and accuracy are reached. 

It appears therefore natural to avoid a partial wave representation completely and work directly with vector vari- 
ables. This is common practice in bound state calculations of few-nucleon systems based on variational [lol | and 
Green's function Monte Carlo (GFMC) methods [H IH [H 111 , which are carried out in configuration space. 

Our aim is to work directly with vector variables in the Faddeev scheme in momentum space. In earlier work |l5l 
we showed that the bound state Faddeev equation has a rather transparent structure when formulated with vector 
variables compared to the coupled set of two-dimensional integral equations obtained in a partial wave decomposed 
form. Based on Malfliet-Tjon type interactions for two- as well as three-body forces, it was demonstrated that the 
numerical solution of the bound state equation using vector variables is straightforward and numerically very accurate. 
As far as three nucleon scattering is concerned, the neutron-deuteron break-up process has been successfully studied 
up to 500 MeV projectile energy based on the first order term of the Faddeev equation using realistic nucleon-nucleon 
forces 0. 

In this article we want to show that the full solution of the three-body scattering equation can be obtained in a 
straightforward manner when employing vector variables, i.e. magnitudes of momenta and angles between the mo- 
mentum vectors. As a simplification we neglect spin and iso-spin degrees of freedom and treat three-boson scattering. 
The interactions employed are of Yukawa type, and no separable approximations are involved. The Faddeev equation 
for three identical bosons is solved exactly as function of momentum vectors below and above the three-body break-up. 

This article is organized as follows. Section II reviews the Faddeev equation for three-body scattering in momentum 
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space and discusses our choice of momentum and angle variables for the unknown amplitude in the equation and its 
kernel. In Section III we derive the amplitudes and cross sections for elastic scattering and break-up processes. In 
addition we relate both via the optical theorem. In Section IV we discuss the numerical methods necessary for solving 
the Faddeev equation, especially our treatment of the singularities in the free three-body propagator. In addition, our 
numerical tests for the solution are shown and discussed. In Section V we present calculations for elastic scattering 
and break-up processes in the intermediate energy regime from 0.2 to 1 GeV. Our focus here is on the study of the 
importance of rescattcring terms as function of the projectile energy and the reaction considered. We conclude in 
Section VI. 



II. FADDEEV EQUATIONS FOR THREE BOSONS IN THE CONTINUUM 

Various presentations of three-body scattering in the Faddeev scheme are presented in the literature 0, 0, E| ■ We 
solve the Faddeev equation for three identical particles in the form 

T\<j>)^tP\(j>)+tPGoT\(j>). (2.1) 

The driving term of this integral equation consists of the two-body t-matrix t, the sum P of a cyclic and anticyclic 
permutation of three particles, and the initial state |</)) = \tpdCio), composed of a two-body bound state and the 
momentum eigenstate of the projectile particle. The kernel of of Eq. (|2.1() contains the free three-body propagator. 
Go = [E — Ho + ie)~^, where E is the total energy in the center of mass frame. 
The operator T determines both, the full break-up amplitude 

C/o = (1 + P)T, (2.2) 

and the amplitude for elastic scattering 

U = PGo ^ + PT. (2.3) 

In this paper we focus on three identical bosons and use a momentum space representation. For solving Eq. (|2.1|) , we 
introduce standard Jacobi momenta p, the relative momentum in the subsystem, and q, the relative momentum of 
the spectator to the subsystem. The momentum states are normalized according to (p'q'|pq) = S^{p' — p)i5'^(q' — q). 
Projecting Eq. (|2.1() on to Jacobi momenta leads to [T^ 

^ 2 A^n J E-^^{q^ + q"^ + q-q")+ie ^ ' 

Here ts(p', p) = t(p, p') + t{—'P', p) is the symmetrized t matrix and E is the total energy in the center of mass (cm.) 
system 

E^Ed + -^ql =Ed + ^Eiab- (2.5) 

We assume that the underlying force is a two-body force, generating t via a two-body Lippmann-Schwinger equation 
and supporting one bound state with energy E^. Thus, ts{z) has a pole at z = Ed- Since the transition operator T 
in Eq. (|2.4|) is needed for all values of q, one encounters this pole of tg . Extracting the residue explicitly by defining 

, 4(p',p,z) 

ts(p,p,2:) = — (2.6) 

z — i^d 

and similarly for T, Eq. (|2.4|l can be rewritten as 

(pq|f |(/?dqo) = (/3d ^q+ iqo^ is ^p, ^q + qo,-B- 

i(p,lq + q",i^-^g^) (q+|q",qlT|^.qo^ 
^ E-^{q^+q"^+Ci-c^') + ie E-^q"^-Ed + ie' 
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This expression is the starting point for our numerical calculation of the transition amplitude without employing an 
angular momentum decomposition. 

The first important step for an explicit calculation is the selection of independent variables. Since we ignore spin 
and iso-spin dependencies, the matrix element {pq\T\tp(icio) is a scalar function of the variables p and q for a given 
projectile momentum Qq. Thus one needs 5 variables to uniquely specify the geometry of the three vectors p, q and 
qo, which are shown in Fig. 1. Having in mind that with three vectors one can span 2 planes, i.e. the p-qo-plane and 
the q-qo-plane, a natural choice of independent variables is 

P = IpI, q = |q|, Xp^p- qo, Xq = q - qo, (qo x q) • (qo x p). (2.8) 

The last variable, is the angle between the two normal vectors of the p-qo-plane and the q-qo-plane, which are 
explicitly given by 

, — X qo X p 
(qo X p) 



v/i-(qo-p)2' 

(qo-^q) = (2.9) 
Vi - (qo ■ q)^ 

It should be pointed out, that the angle between the vectors q and p, ypq = p • q, is not an independent variable. In 
fact, if Xp and Xq are given, the domain of ypq is bound by 



1 - 2:2^1 ~xl< ypq < XpXq + Jl - xl-Jl - x\, (2.10) 



thus not covering the entire interval [-1,1]- Using the explicit representation of the normal vectors and standard cross 
product identities, we arrive at the following relation between x'^q and ypq. 



p ■ q - (qo ■ p)(qo ■ q) 



- (qo • P) Vi - (qo • q)' 



(2.11) 



Upq ~ XpXq + Y 1 -^py^ •'^'q ■^'pq' (2-12) 

For the special case where qo is parallel to the z-axis ((/o-system) one can write 



ypq = XpXq + -y/l - 3;2y^l - X^ COS(ppq, (2.13) 

where the ippq is the difference of the azimuthal angles of p and q. However, the variable cos ippq, which was used 
erroneously in |l9l | as third angular variable, is not rotationally invariant. 

With the independent variables listed in Eq. H2.8|l the matrix element of T is given as 

(pq|r|^dqo) = f (p, Xp, xll, Xq, q; qo). (2.14) 

Furthermore, ts{p' ,p, z) is also a scalar function, and thus can be written in the form 

i(p',p,z) =i(p',p,p'-p,z). (2.15) 

The most intricate dependence appears under the integral in Eq. H2.7|l for the third angular variable of the T- 
amplitude. According to Eq. 12.11|l it is given as 

_ (q^q^O-q"-qo-(q^q'Oqo-q^^ . 

1- qo-(q+H") - (qo • q")' 
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In view of the break-up singularities of the first denominator in Eq. (|2.7() it is mandatory to choose the coordinate 
system for the q" -integration such that the z-axis points paraUel to the vector q. Then one obtains for Eq. (|2.1()|) 



qx" + ^q" 



g+j?" V q 



= q" • qo = x"Xq + a/I - X""^. 1 - Xl COs{ip" - LPqg) 



where 



1 X QXg + \q"x„ii 

Xq^y, ^ (q + -q") • qo = ' ' " =. (2.18) 

Here ipq,^ is the azimuthal angle of qo in the coordinate system chosen for the q" -integration. These considerations 
lead to the explicit representation for the transition amplitude T 



f{p,Xp,xf> Xg,q;qo) = (fid i \ + ^Qo + qqoXg] (2.19) 



X U I P, \l -y + % + moxq, =; E~—q 

\/ig'+go+ggoa;g *™ 

f + l <.2ir 1 



E- Uq^ + qq"x" + q"^)+ie 



-( I 1^2, 111 \ II II 2 gyP9 "'^ 1 VPI" TP 3 2 

X i» P. A/ 7g + g + 2;", . ;-E - —q 

4 . /I ^2 + qii2 + ^q.^,,, 4m 



T 



Iq2 + i ia + qqii^ii^ „ , ^^'^-p^'^;:^ ,x,..,q";q, 

^-4lrg"'-^^' + *^ 

where in addition to Eq. H2.8|l and the related variables of Eq. H2.12|) and Eq. H2.18|l the following variables occur 



g" = |q"l, 

x" = q-q", 



Vpq" = P • q = Vpq- 



c" + v/l-x"yi-y2^cos(v.p-^"). (2.20) 



Like ifq^ in Eq. H2.18|l . the angle ipp in Eq. (|2.2U|) is the azimuthal angle of p in the g-system (i.e. the system where 
the z-axis is parallel to q). It remains to relate the angles Lpp and (^gj, to the three angular variables Xp, Xq, and . 
As it is shown in Appendix 1X1 due to the (y9"-integration, only the knowledge of cos(iy9p — ^Pq^) is required. Like cos ifpq 
in Eq. (|2.13() is equal to in the gg-system, so is cos((/3p — (pq^^) equal to x^^^p in the g-system. Thus 

I . q go -p - (q-qo)(q-p) xp - Xgy^q 

cos{(pp -^ipq,) = xl = ^ = = ^ ^=. (2.21) 

\/l-(q-qo)Vl-(P-q)' Jl^x^„ Jl-yL 



Because of that difference ifp — ^qo)i one can choose (pq^ arbitrarily, e.g. zero. Furthermore, cos (pp and sin ipp required 
in Eq. (|2.20() are given in terms of cos((/9p — (fqa), as is shown in Appendix fXl This completes the definition of all 
relevant variables in Eq. (|2.19f) . 
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III. AMPLITUDES AND CROSS SECTIONS FOR ELASTIC SCATTERING AND BREAK-UP 

PROCESSES 

The amplitude for elastic scattering is obtained by calculating the matrix element of the operator U given in 
Eq. lESJl as. 

(c[ipd\U\qofd) = 2¥3<jQq + qoj f^' - ^(g^ + q • qo + ^o) ] V'd + ^^o 



r ^ n \ (q + iq",q"|r|qo^d; 

The amplitude for the full break-up process according to Eq. ()2.2|l . is given by 

fpq|f|qo</Jd) (-5P+|q,-p-iq|r|qo¥'d) (-5P - |q, +P - |q|r|qo</'d 



(pq|[/o|qo<(5d) 



(3.2) 



The equation for the elastic operator U follows from Eqs. I|2.1|l and 12.3|l . It is given as 

f/|(/)) = PGo ^ |(/.) + PtGoU\q}) (3.3) 
Straightforward and well known steps |^ based on this equation lead to the unitarity relation 

(0|C/|0r - = J d'q{<j>,\U\cl,r27Ti6{E~Ed^^q'){c^,\U\<j>) 

^ d^pd'q{MUom*2^^S{E-^-^q^){^o\Uo\^)- (3.4) 



3 / m 4m 



We want to point out that there is a misprint in Eq. (202) of Rcf. the factor 1/3 is missing. 

Using the variables defined in the previous section, and having in mind that for elastic scattering |q| = |qo|, the 
amplitude for elastic scattering according to Eq. H3.1|l can be expressed as 



{q^d\U\qoVd) = U{qo,Xg)^2^l(^qo^l^+7gj (^E - ^{2 + Xq) ] (3.5) 
+ 2 / dq"q"^ / dx" / d^" 



J-1 Jo E^^q"^-Ed + ie 



X 'Pd\\l ^ql + + qoq"yqq" 



X T 



2 n1 , // q0Xq + ^q"yq„q" ^ ql + \q"^ +qm" V,," '''' „ 

% + 1^ + <iaq"yqq", , I 7 > %o«". 9 ; 90 



y y q5 + ji + loQ Vqq 



with 



Vqq" = q-q" 
Vqoq" = qo • q" = 



qgXq + \q"yqaq" 

ql + \q"^' + qaq"yqq" 



(3.6) 



At this point, the choice of a specific coordinate system for the q"-integration is still open. The angular variable 
= q • qo represents the scattering angle. If the z-axis is chosen parallel to qo, the angles are 

Vqq" = XqX" + ^ 1 - X^^/l - x"^ COa{(pq - (fi"). 
Vqoq" = (3-7) 
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While having the z-axis parallel to qo may be the intuitive choice to describe the scattering with a given beam 
direction, wc can in principle also choose the z-axis parallel to q. In that case the angels in Eq. I|3.6|) are given by 



Uqq" X 



Vqoq" 



" + ^l~xl^l-x"^ cos{^g„ ~ ^"). (3.8) 



The elastic cross section depends on the angle between the vectors qo and q, but not on the choice of z-axis. We use 
the possibility of calculating the matrix elements of U in the two different coordinate systems to test the quality of 
our numerical calculations. 

The differential elastic cross section in the cm. frame is given by 

2 



— (27^n[/(gox,)|^ (3.9) 



and the corresponding total elastic cross section is 

Oei= J ^ (^^y {2nf j''\x\\U{qo,xr. (3.10) 

The full break-up amplitude is given in Eq. (|3.2|) . On the energy shell p and q are constrained by + jq^ 
As function of all five variables and the projectile momentum it reads 

, qg . _ T{p,Xp,x'^°g,Xq,q,qo) T(p2, 

uo{P,Xp,Xp„,Xq,q,qQ) = — n - — I -:; 3 5 -:; I -:; 7^ ■ I'^'-'^-'^J 



E - 32-^2 _ ^ E- -^ql -Ed E- -j2-g2 _ £; 

Here the variables are defined as 



Upq — XpXq + a/1 Xp-ul Xg 



I 1 3 , 1 / „ 9 „ ^ 
P2 = \- + ^(l\ = -dp'^ + -q^ -3pqyp 



92 = I - P - ^q| = Y + + PqVpq 



1 3,1/^9 



Pa = 1^2^^ i'^' " 2 V^^ ^ 4*^^ ^ ^PlVpq 



93 = I + p - ^q| = Y + ^9^ - pqvpq 

1 3 
yP-^p ~l" iQ-^q 

= P2 ■ qo = — — 

P2 

T:q2 = q2 • qo = — 

92 

1 3 
^ * ^2P'^P '^Q'^q 

^P3 — P3 • qo — - 



qs • qo 



P3 



93 

40 = (qo'^P2) • (qo'^qs) = ^ 



(qo X P3) • (qo X qs) = (3.12) 

l-xl^Jl-xl 



P3 y 93 

The five-fold differential break-up cross section is given in the cm. frame 



d^abr (27r)'^m2 2 / ^ - 2 irr / 90 m2 /oiqa 

- -qi mE--q-' \Uo{p,Xp,x'i" Xg,q,qo)\ (3.13) 



dflpdflqdq 3qQ 
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It is convenient to calculate the total break-up cross section in the cm. frame, since there are no kinematic 
restrictions on the relative angles. For the explicit calculation we can make different choices of the z-axis, e.g. it can 
be parallel to the direction qo of the projectile, or parallel to either one of the Jacobi vectors q and p. The different 
choices will obviously result in different angular integrations. For completeness we give all three choices here. This 
will be used as a highly non-trivial test of the numerical results, as will be demonstrated in the next Section. If the 
z-axis is parallel to qo we have 

27r 



-1 

with 

Xp > Xp , Xq > Xq , "^pq 



d^pdVLq ^2ti I dxl I dx'^ J difi'^g (3.14) 

cosippq. (3.15) 



If the z-axis is parallel to q, the angular integration becomes 



+ 1 r+l I-2TV 

dVLpdQ.q = 2iT j dxg J dy'pq I d(fpq^ (3.16) 



with 



XqVpq + Jl - X'^jl - yfq COS ^^'^^ , 



1 - [x'^Vpq + ^1 - ^1 - 2/;'2 COS^'^q^^ 

Finally, if the z-axis is parallel to p, the angular integration is 

f + l /•27r 



with 



Ojp 
Xn 



Xpy'pq + ^^-^f\j^-yfq COS '^"qqo ' 

y'U - {^'ly'pq + y^i -xf^\-yf^ cos (^'^'^, 



x" 

J^p 



1 - ( + \/i - - y'm ^^^vqqo \/i - 



(3.17) 



J dnpdQq = ^^J^ dx'; J ^ dy';^ ^ (3.18) 



(3.19) 



T{p,q) ^ / d%dnq\{cj)o\U^m^ , (3.20) 



Let us define a function !F{p^ q) as 



where the angle integrations over the break-up amplitude is carried out. This function should be independent of 
the coordinate system in which the angle integrations are performed. We use this property to check our numerical 
calculations. This is a non-trivial test of our calculation, since especially at higher energies the transition amplitude 
T develops stronger angle dependencies, which challenge the accuracy of the multi-dimensional interpolations. 
The angle integrated break-up cross section is given as 

^ - ^^^'^^Z"^^- ^ I ' (3.21) 

and the total break-up cross section reads 




l(27r)4m2 /-V^ 



cJbr^T^ ' / dqq\l7nE--q^ T\xlmE~-y,q]. (3.22) 
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Using now the unitarity relation from Eq. H3.4() the optical theorem gives 

4m(27r)'^ 

o-tot = (Jei + cTfcr 5 Imt/ (go, 1) • (3.23) 

For later use we also mention the semi-exclusive cross section, where only one particle is detected in the break-up 
process 



da 



dflqdq 



{27T)^—pq^ J dp\Uoip,Xp,x;';j,Xg,q,qo)\'' (3.24) 



IV. NUMERICAL METHODS 



The fully off-shell two-body t-matrix t(p', p, z) is solved directly from the Lippmann-Schwinger equation as function 
of its vector variables for the off-shell energies E— -£^q^ as required by Eq. H2.4|l . The Faddeev equation is iterated, 



generating the multiple scattering series, which is then summed by the Fade method |2lll22l |. We use it in the form 
of a continued fraction expansion as layed out in Ref. [l^ . 

The first integration to be performed in solving Eq. (|2.19|l by iteration is the integration over the azimuthal angle 
ip" . This leads to a function of variables, q" and x" and requires interpolation in the second and third arguments of 
ts and the first four arguments of T. Our spline interpolations are based on the cubic Hermite splines given in 231 ■ 

Let F{q",x") be the resulting function in each step of the iteration. Clearly, it depends in addition on the fixed 
variables p, q, Xp, Xq, and x'^°q, which are omitted for clarity. Then the next task is performing the remaining two 
singular integrations, 

dn"a"^ dx" F{q"^x") 

J -I {E-^^{q^ + g"2 + qq"x") + le) {E - -^q"^ -Ea + ie) 

If the cm. energy E is below the three-body break-up threshold, only the second denominator is singular, and the 
simple pole can be treated by standard subtraction methods. 

The intricate problem arises above the three-body break-up threshold, when in addition the first denominator can 
vanish. It has the form 



: = (4 2) 

E- ^^{q^+ qq"x" + g"2) + ie Xq - x" + ie ^ ' ' 

with 

mE-q^-q"^ 

xo = . 4.3 

qq" 

For |xo| < 1 a so-called moving singularity arises in the g"-a:"-integration, since depends on q. The direct 
treatment of those moving singularities using real variables has been discussed in the literature 0. We briefly review 
the appearance of these singularities in form of logarithms, since we introduce a new quasi-analytic integration based 
on Spline functions. The condition Ixol = 1 leads to the pole positions 



q" =±^±^mE-\q^ (4.4) 
and one arrives at the well known shape in the (/-(/"-plane for |a;o| < 1 shown in Fig. 2. This region is bounded by 



q+ = l + jQl-\q^ (4.5) 



and 
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where Qq — V mE. Apparently there is no singularity if g > Qmax = y We distinguish four cases, q = 0, 

< q < qmax, q = qmax, and q > qmax- The case g = reduces to a simple subtraction and will not be discussed. 
For < g < qmax we only consider the part of the (/"-integration which contains the moving singularities. It has the 
schematic form 



dq" 



+1 



dx ' 



f{q".x") 

xq — x" + is 



(4.7) 



The first step is to perform a subtraction of the pole, which we carry out in the entire square < g, g" < qmax by 
defining 



We obtain 



fiq",Xo): \xo\<l 
/(g",^): |xo|>l 



dq" 



1 



rfg'7('7",^o)ln 





Xo - 


-x" 


1 


+ Xq 


. f 


1 






- Xo 


Jo 



dq"e{i-\xo\)nq",xo). 



(4.8) 



(4.9) 



Now we define g_ 



In 



1 


+ Xo 


-{ 


1 


- Xo 





Qq — |g2 and obtain 

q 



ln|g" + |g_||-ln|g"+g. 



ln|g"-|g_||+ln|g"-g^ 



This leads to the well separated part of the integral which contains the logarithmic singularity 

'"'""^ dq"f{q",Xo) In |g" - |g-|| + In |g" - g, 

V 19-1 



(4.10) 



(4.11) 



It is here that we introduce the new technique that relies on cubic spline integration. 

We divide the range of integration [0, qmax] into intervals bounded by a set of grid points g^. The set of grid points is 
supposed to be dense enough to interpolate the function f{q", xq) = f{q") sufficiently well by cubic Hermitean splines 
[231. In Ref. a detailed presentation of these spline functions is given specific for our use. For the convenience of 
the reader we now switch to the notation of |23l | and denote the end points of the i-th interval by Xi and X2 and the 
two adjacent grid points to the left and right of the i-th interval by xo and X3, respectively. Then as detailed in 
the interpolating function in the i-th interval (replacing f{q") = fix)) can be written as 



fiix) = ^S'j(a;)/(a;j), 



where the modified spline functions are 

X2- Xi 



with 



So{x) 
Si{x) 
S2{x) 
Ssix) 



-03 (x)- 



1 



(/)i(a;) + 03 

(j)2{x) + 03 
(j)i{x) 



Xi ~ Xo X2 — Xo 

X2 — Xl Xi — Xo 



Xi — Xo X2 — Xl J X2 — Xo 



4>iix) 



X3 ~ X2 



X2 — Xl X3 — Xl 



Xl - Xo 



1 



X2 - X1X2- Xo 
X2~ Xl 1 



+ 4>4.ix) 



X3 ~~ X2 X2- Xl 



X2 — Xl X3 — X2 J X3 — Xl 



X3 - X2 X3~ Xl 



Mx) 



ix2 - X)'^ 
(X2 - Xi)^ 



[{X2 - Xl) + 2{x - Xl) 



(4.12) 



(4.13) 
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<l>2{x) = 7 [(x2 - xi) + 2(a;2 - x)] , 

[X2 - Xi)-^ 



2 



{X - Xi){x2 - X ) 

{X2 - XiY 
{x — Xi)^{x — X2) 
{X2 - XlY 



M^) = ' rJ' ":,2 ■ (4-14) 



Therefore, in view of Eqs. (|4.11() - (|4.14|l . the following integrals occur for i = 1 • • • 4 

4>i — I (f)ihi\x — q\dx, 



/ 02 In |x - qldx, 

J Xi 

Xi + l 

03 In |a; — q\dx, 

Xi+l 

04 In |x - 5|da;. (4-15) 



Xi 



with q — g+. Consequently the five different cases q < Xi < x^+i, q = Xi < x^+i, Xi < q < Xi+i, Xi < q = Xi+i, 
and Xi < Xi+l < q occur. Since the functions 0i(x) are cubic polynomials, the integrals in Eq. H4.15|l can be performed 
analytically. We leave the explicit calculation to interested practitioner and refer to for a detailed presentation. 
According to our experience that manner to integrate the moving logarithmic singularities is a very good alternative 
to the more common subtraction method 01- 

Finally, for q = qmax we also apply the subtraction over the extended region < q" < qmax- In that case 
q+ = g_ = and when q" = then xq = — 1. Analogous steps lead to that part of the integral, which contains 
the logarithmic singularity, 

"""^^'^ dq"f{q",-l)\n\q"-^\, (4.16) 

which is again evaluated by spline based integration. 

In order to test the correctness as well as the accuracy of our calculations we carried out a variety of numerical 
tests. Unfortunately we could not compare to work by other groups, since to the best of our knowledge no comparable 
work at higher energies exists. 

Apart from the projectile momentum qo, the amplitude T of Eq. (|2.19|) depends on five variables p, Xp, x'^''^^Xq,aiid 
q. In addition, there are the integration variables q", x" and ip" . All calculations listed are based on the Malfliet- 
Tjon type potential which is explicitly given in the next section. The fully off-shell two-body i-matrix, t(j)' ,p,x,e), 
is obtained for each fixed energy on a symmetric momentum grid with 60 p {p') points and 40 x points. Since 
the momentum region which contributes to a solution of the two-body t-matrix is quite different from the region 
of importance in a three-body calculation, we map our solution for tg onto the momentum grid relevant for the 
three-body transition amplitude. This is done by applying the Lippmann-Schwinger equation repeatedly. The t- 
matrix ts{p' ,p, x,e) is obtained at energies e — E — ■^Q^, exactly at the q values needed in the three-body transition 
amplitude of Eq. 

In carrying out our calculations, it turns out that there are essentially two separate issues governing the quality 
of the results. The first is the angle dependence of the transition amplitude of Eq. (|2.19|l . It is to be expected that 
the angle dependence is weak at low energies and increases with higher energies, reflecting the need to include more 
and more partial waves at higher energies in a partial wave based calculations. As example we list in Table n the 
elastic and break-up total cross sections together with the total cross section extracted from the imaginary part of U 
in forward direction, Eq. H3.23|l . At 0.01 GeV 12 points for all angles are clearly sufficient, whereas at 0.1 GeV this 
is not so. Tabled lists the elastic, break-up and total cross section as function of the angle variables, and we see that 
one needs at least 16 points for all angles. At 0.5 GeV we find that the biggest angular dependence occurs in Xg and 
x" , and the least dependence in the azimuthal angle and the angle x^° , and take this into account in our choice of 
angle points. 

The other issue is the quality of the calculation in the singular regime, i.e. in the integration region bounded by 
qmax in Fig. 13 We divide the integration grid for g" into the intervals {0,qmax) U {qmax.q), where q^ax = 
and q =20 fm^^. The inteval boundaries and qmax are handled explicitly. As the energy increases, qmax increases, 
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and we need to take this into account by changing the distribution of the g" points as function of energy within the 
different (j"-intervals, i.e. put more points into the interval (0, Qmax) and less into {qmaxj q)- From the number of points 
in {0,qmax) one can define an average point distance = (7^02; /(number of points m{0, qmax j) in this interval. In 
Fig.|21we show the dependence of the calculation on Aq by using the percent error Sgpt = {(Jopt — o'ei — o'brl/f^opt x 100 
in the fulfillment of the optical theorem as quality measure. At 0.01 GeV it is quite easy to make the average point 
distance very small in the interval (0, qmax), since qmax is only 75 MeV. The top panel of Fig.|31shows that the percent 
error 5 opt drops linearly below 0.1% and flattens out at Ag= 3.5 McV where most likely errors in the interpolation 
start to play a role. At projectile energy 0.1 GcV qmax is already 284 MeV, and is naturally much larger with 
a reasonable number of q" points. The dependence of the 5opt for 0.1 GcV on Ag is shown in the middle panel of 
Fig. 01 for two different cases. An angular grid size of 12 points is indicated by the open squares, one of 24 points by 
crosses. A comparison of the calculations shows at Ag = 15 MeV the calculation with 12 angular points can not be 
improved any further, the calculations start to oscillate for smaller Ag, accidentally giving a very good agreement at 
Ag = 13 MeV. Increasing the number of angle points to 24 shows a further linear decrease in the error (cross symbols) 
into the 1% region of 5opt at a Ag = 9 MeV. This is consistent with the findings shown in the top panel. We continue 
to study the dependence of Ag at 0.5 GeV, where qmax = 644 MeV. Here we immediately use 24 angle points as 
suggested from Tabic A total of 30 integration points in {0,qmax) leads to Ag = 22 MeV and Sopt ~ 10%, which 
is consistent with values in the middle panel. The 10% error is also consistent with the value for the total break-up 
cross section in the last row of Table which indicates that abr is not yet converged. From the systematics at the 
different energies shown in Fig. |2| we can extrapolate on the Ag needed to reduce the error in the optical theorem. 
Due to computer time limitations we have not pushed this any further. 

In Tabled the total elastic, the total breakup as well as the total cross section evaluated according to Eq. (|3.23|) 
via the optical theorem are given and shown as function of various sets of grid points. The momentum grids for p 
and q arc discretized with 49 points each. The integration variable q" plays the same role as q, and is therefore also 
discretized with 49 points distributed over the intervals {0,qmax) U {qmax,q) in an energy dependent way according 
to the insights described above. The values given in last row of each energy correspond to the points at the smallest 



A nontrivial test for the quality accuracy of our calculation is the numerical verification of the optical theorem 
Eq. (|3.23() . Our results are given for selected energies in Table Hll Here we show two sets of cross sections, distinguished 
by the superscripts go and q for the total and clastic cross section. The superscripts indicate that the calculation is 
carried out by choosing the z-axis cither parallel to qo or to q. Performing the calculation with two different choices 
of the z-axis is a nontrivial test for our choice of independent variables as well as for the entire calculation. The 
total break-up cross section is also calculated in a coordinate system in which the z-axis is parallel to p, indicated by 
a^^. The calculations are based on the largest grids given in Tableland show a very good agreement of the results 
obtained in the different coordinate systems. This indicates the numerical rotational invariance of our calculations. 

On top of convergence tests for the Fade summation, we insert the resulting amplitude T again into the integral of 
Eq. H2.19|l . leading to a second amplitude T' . Both amplitudes should be identical within our numerical errors. This 
is checked by evaluating the cross sections again using the second amplitude. We document the results in Table Hill 
for the differential clastic cross section at selected angles. The tabic shows excellent agreement of the two values of 
the cross section 

Finally, another highly nontrivial test of our calculation is the independence of the cross sections from the arbitrary 
angle ipg^ and the sign of sm{ipp ~ fq^). This is documented in Tables HVl and Ivl for the energy £^=3 MeV. In order 
the check the rotational invariance numerically, the calculations are carried out in two different coordinate systems, 
one where qo is parallel the z-axis, and one where q is parallel to z. Both tables show excellent agreement for the 
cross sections, and thus we conclude that our choice of variables is correct. 



Through we neglect spin and iso-spin degrees of freedom and stay in a strictly nonrelativistic framework, we 
nevertheless can provide first qualitative insights for various cross sections in three-body scattering in the intermediate 
energy regime which we define from 200 MeV to 1 GeV projectile energy. The focus of our investigations will be the 
question which orders of rescattering in the two-body i- matrix are needed to come close to the exact result, namely 
the solution of the Faddeev equation. 

As a model two-body interaction wc choose a superposition of two Yukawa interactions of Malffiet-Tjon type 



Ag in Fig.El 



V. 



SCATTERING CALCULATIONS AT INTERMEDIATE ENERGIES 




(5.1) 



The parameters are given in Table IVII and fitted such that the potential supports a two-body bound state, the 
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'deutcron', at -2.23 MeV. As first result we show in Fig.^jthe total cross section together with the total elastic and 
total break-up cross section as function of the laboratory projectile energy. In addition, the total cross section is also 
evaluated via the optical theorem as a test of our numerics. This duplicates the information already given in Table III. 
We see that the optical theorem is quite well fulfilled. The figure shows that at roughly 1 GeV the total elastic and 
total break-up cross section become equal in magnitude in our model. 

Next we show in Fig. 4 the angular distribution in elastic scattering for a set of selected energies. In addition to the 
exact Faddeev result, the cross sections are evaluated in first order in the two-body t-matrix, second order in t, third 
order in t, and 4th order in t, and displayed. First we notice that with increasing energy the cross section in forward 
direction decreases. Furthermore, for all energies shown, the first rescattering (2nd order in t) always increases the 
cross section, and subsequent rescatterings lower it again. As expected, for the lowest energy, 0.2 GeV, rescattering 
terms of higher order are important, and even the 4th order is not yet close to the full result. The same is true for 
0.8 GeV. We notice, that even at 1 GeV two rescattering terms (3rd order in t) are necessary to come into the vicinity 
of the final result. The same is true for 0.5 GeV. 

In view of the standard '^-p' impulse approximation for the optical potential in nucleon- nucleus scattering employed 
at intermediate energies |26j. it is interesting to notice that the first order result in t in our model study is quite 
insufficient. Even at energies larger than 0.5 GeV rescattering corrections up to the 3rd order arc required to come 
close to the exact result for small scattering angles. Therefore it seems to be likely that the first order impulse 
approximation in nucleon- nucleus scattering is insufficient. 

In case of inelastic processes we first regard the semi-exclusive reaction d(N,N') where only one nucleon is detected. 
We choose three different laboratory energies, 200 MeV, 500 MeV, and 1 GeV and show the inclusive cross section as 
a few selected angles for the detected nucleon. The results are shown in Figs. 1^ through [TUl 

At 0.2 GeV the semi-exclusive cross section is given in Fig.|Hlfor the emission angle 24° and in Fig.[7|for the emission 
angle 39°. Both figures show in the upper panel the entire energy range of the emitted particle. Since the cross section 
varies by two orders of magnitude, we display it in a logarithmic scale. In order to better flash out the peak structures, 
the two lower panels show the high and low energies of the emitted particle in a linear scale. Together with the full 
solution of the Faddeev equation (solid line) we display the sums of the lowest orders of the multiple scattering series 
as indicated in the figure. The peak at the highest energy of the emitted particle is the socalled final state interaction 
(FSI) peak, which only develops when rescattering terms are taken into account. This peak is a general feature of 
semi-exclusive scattering and is present for all energies. The next peak is the quasi-free (QFS) peak, and one sees that 
at both angles one needs at least rescattering of 4th order to come close to the full result. However, in contrast to the 
smaller angle, at the larger angle, 39° in Fig.[7| the first order result for the larger energies is surprisingly close to the 
full solution, though the multiple scattering series is by no means converged, as the following higher orders indicate. 
We also observe that the QFS peak moves to lower energies of the emitted particle with increasing emission angle. 
At both angles the very low energies of the emitted particle exhibit a strong peak in first order, which is considerably 
lowered by the first rescattering. Here the calculation up to 3rd order in the multiple scattering series seems already 
sufficient. 

For 0.5 GeV incident energy the semi-exclusive cross section is given in Fig. |Slfor the emission angle 24° and in 
Fig. El for the emission angle 36°. We again see three peaks along along the energy axis of the detected nucleon, 
the FSI and QFS peaks as well as the peak at the extreme low energy of the emitted particle. Again we see that 
the results based on first and second order in t alone are quite insufficient and higher order rescatterings can not be 
neglected. It is also interesting to observe that at 24°, Fig.|Sl the third and higher order rescattering terms shift the 
peak to higher energies, whereas at the larger angle of 36° the peak positions of the various orders coincide more 
or less and agree with the peak position of the full calculation. Again, for the peak for the very low energies of the 
emitted particle, the 3rd order calculation agrees already quite well with the full result. 

At 1 GeV the situation is similar. As examples we have selected two angles, 18° displayed in Fig.^|and 30° displayed 
in Fig. ^] For the small emission angle the second rescattering shifts the QFS peak towards higher energies, at the 
larger angle this is not the case. Our studies indicate that this is a general phenomemon occurring at all energies 
under consideration. There exists a critical maximum energy E^""^ of the emitted particle, corresponding to a specific 
emission angle, at which such a shift in the QFS peak through higher order rescattering terms can occur. At 0.5 GeV 
projectile energy this maximum energy is 0.44 GeV, at 1 GcV it is 0.88 GeV. If one considers the ratio -g — , then 
one finds for both cases that if this ratio is larger than 0.8, the QSF peak is shifted by higher order rescattering terms. 
This could be interpreted as an interference between the QFS and the FSI mechanisms. If this ratio is smaller than 
0.8, then the FSI peak is small and higher orders in the multiple scattering series do not change the position of the 
QFS peak. In addition, it seems that at the larger angle (Fig. Ill|l the multiple scattering series converges a little 
faster with respect to the higher orders compared to the smaller angle fFig. I10|l . The final result for the peak at the 
very low energy of the emitted particle is as before reached with two rescattering contributions. It is remarkable that 
for the energies between about 200 and 500 MeV of the emitted particle the first rescattering contributes almost an 
order of magnitude to the cross section. 



13 



We can make a first contact to calculations based on realistic nucleon-nucleon (NN) forces. In Ref. [T3 the semi- 
exclusive process d{p, n) has been determined in first order in t based on the NN potentials AV18 [23| and Bonn-B 
[2^ . In the upper panel of Fig. El we compare our first order calculation at projectile laboratory energy 495 MeV 
and 18° emission angle with the first order calculations from Ref. ^3 based on the two realistic potentials. The 
position of the peak is only determined through kinematics, thus the peak position coincides for all three calculations. 
Though our model calculation refers to bosons and the potential contains only the crude features of a central short 
range repulsion and intermediate range attraction, the magnitudes of the cross sections differ only by roughly 20%. 
In the lower panel wc show the contributions of the first orders of the multiple scattering scries successively summed 
together with the exact solution of the Faddeev equation for our model. At this angle and energy the contribution 
of the first rescattering (2nd-order in the multiple scattering series) is quite weak, the contributions of the next two 
orders are large and lower the size of the peak. At the very high energies of the emitted particle, the 4th order in 
the multiple scattering series is still not yet close to the exact result. Therefore we conjecture that at this energy 
calculations with realistic forces will also require higher order rescattering contributions. 

Finally we comment on a recently measured and analyzed reaction pd — > {pp)n at high momentum transfer pfll l3nl| . 
In this experiment the break-up configuration has been chosen such that the neutron is ejected at extreme backward 
angles, and the two protons at extreme forward angles. The measurement was carried out at GcV laboratory energies. 
The experimental data [23 have been analyzed in Ref. [s^] using first and second order processes in the NN t-matrix 
including a A-isobar mechanism. Within our nonrelativistic toy model for three bosons we are of course unable to 
analyze the data. However, within our model we can give a clear answer whether higher order rescattering processes 
are essential in this reaction. In |29ll30l | the data are integrated over a small interval of the relative pp energy between 
and 3 MeV, and averaged over the neutron cm. angle in the interval between 172° and 180°. In our qualitative 
study we fix the cm. angle of one particle at 180°, but integrate over the relative energy Epp = p^ jm of the two 
other particles between and 3 MeV. Thus, wc evaluate the cross section 



da (2\^m^ f^fint:^,^ ^ 

— = (27rf (^-j — dpp'qj dp\Uo(p,Xp,x%,Xq = -l,q,qo?. (5.2) 



Since we choose the z-axis parallel to qo, and q is antiparallel to qo, the ipp dependence is directly given by = cos Lpp. 
Our calculations are carried out for projectile laboratory energies between 0.2 and 1 GeV, and are displayed in Fig. 1131 
Here we compare different low orders in the two-body f-matrix with the full solution of the Faddeev equation. We 
notice that all our calculations exhibit a smooth fall-off as function of the projectile energy. This behavior is present 
in the data of Ref. l29|. None of our calculations shows a dip structure around 0.7 GeV as indicated for some of the 
calculations in Ref. |30j. The reason ma y be that our calculation is carried out in three dimension, i.e. all partial waves 
are included exactly, where as in Ref. [331 only the lowest partial waves are considered. At low projectile energies 
rescattering terms of higher order still give considerable contributions to the cross section. At 1 GeV the first order 
calculation is an order of magnitude smaller than the result of the full calculation. It is interesting to notice that at 
1 GeV the contribution from the first rescattering is relatively small, and one needs to go to the 3rd order in t to come 
close to the full result for this particular break-up configuration. 

In addition to the specific break-up configuration described above, a measurement of the extreme backward scatter- 
ing elastic pd cross section has been investigated in [30| . Instead of a forward scattered pp pair with very small relative 
energy, one now has a forward going deuteron. This situation corresponds to elastic scattering from a deuteron at 
backward angle. To investigate the influence of rescattering for this reaction we plot in Fig. E| the backward angle 
of the elastic cross section at energies from 0.2 to 1 GeV, and compare the result of the full Faddeev calculation with 
calculations based on low orders in the multiple scattering series. The figure shows that the first order calculation 
is insufficient over the entire energy regime considered, except of course for the crossing point around 0.5 GeV. The 
first rescattering contribution (2nd order calculation), though close at 0.2 GeV, is insufficient below roughly 0.9 GeV. 
The figure also shows that at about 0.9 GeV the relative magnitude of contributions from the second and higher 
rescattering become small. Thus, we conclude that at 1 GeV one needs at least one rescattering to be in the vicinity 
of the full result for the elastic cross section at the backward angle. 



VI. SUMMARY AND CONCLUSIONS 



In this study we perform fully converged Faddeev calculations for three identical bosons interacting by non-separable 
forces in the intermediate energy range between about 0.2 and 1.0 GeV. To the best of our knowledge these are the 
first calculations of this kind. The key point is to neglect the partial wave decomposition generally used at low energies 
and to work directly with momentum vectors. Thus all partial waves are exactly included. Important is the suitable 
choice of variables. Besides the two magnitudes of the two relative Jacobi momenta p and q we choose the angles 
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between the vectors p and qo and between q and qo, where qo is the projectile momentum. The fifth variable is 
the angle between the two planes spanned by p, qo and q, qo- In the technical piece of the work we introduce for 
the first time a spline based integration of the moving logarithmic singularities, which is a very valuable alternative 
to procedures used so far. The numerical results are converging as documented in Section III. In Section V we show 
elastic and inelastic (break-up) cross sections in the above mentioned intermediate energy range. We focus on the 
question how many orders of rescattering beyond the often used first order calculation in the two-body t-matrix are 
needed to come close to the full Faddccv result. We find that in nearly all cases studied processes of at least 2nd 
and 3rd order rescattering are required. Whether this will be also required in performing calculations with realistic 
dynamical inputs has to be seen in the future. 

In one case we can make first contact to a result based on the NN forces AV18 and Bonn B, which are considered 
to be realistic in the sense that they describe all NN data below 350 MeV extremely well. This was the semi-exclusive 
cross section at 495 MeV evaluated in first order in the NN t-matrix. Of course, at that energy AV18 and Bonn-B 
are at the upper limit of their applicability. Despite our simple two-body model force, a superposition of two Yukawa 
interactions, one attractive the other repulsive, our results turns out to be within about 20% to the calculation based 
on the realistic models. This shows that our investigations might allow some conclusions about results based on 
present and future models with more dynamical inputs. 

As a first example for considering data in the light of our toy model we study the extreme backward elastic dN 
scattering over the energy range from 0.2 to 1.0 GeV. We find that first order results in the two-body i-matrix are 
totally insufficient and only around 1 GeV the first order rescattering comes close to the full result. Parallel to those 
data in elastic scattering in |2^ the complete break up process d(p,n)pp has also been investigated. Here the 
neutron was ejected antiparallel to the beam direction and the two protons at extreme forward angles with a very 
small relative energy. Again we study the significance of rescattering processes, and find that for this particular 
break-up configuration two rcscatterings arc necessary to get close to the result of the full Faddccv calculation. 

In conclusion we can say that the three-body Faddeev equations can be safely solved at intermediate energies using 
directly momentum vectors. Calculations based on partial wave decomposition would be hardly feasible at these 
energies. 

Further studies scanning the complete three-body phase space for the total break-up are underway. This may 
be important in order to shed light on previous theoretical analysis of p(d,ppn) reactions which relied on low order 
reaction mechanisms. 

Based on our current experience it appears that if low order rescattering processes will turn out to be sufficient 
for certain phase space regions, realistic calculations including spin and isospin will be feasible, even including three- 
body forces. A first step evaluating the d(p,n)pp break-up cross section in first order with a currently used two-pion 
exchange three-nucleon force model is already under way [sij . What is badly needed now are realistic models for 
nuclear forces in the intermediate energy regime we study. This paper allows us to conclude that it will be feasible to 
extend the calculations to realistic dynamics. 

In Ref. 01 the effect of relativistic kinematics was studied in the (p,n) charge exchange reaction on deuterium 
between 0.1 and 0.5 GcV in a first order Faddeev calculation. This work concluded that the effects due to relativistic 
kinematics arc quite visible at 0.5 GeV, specifically in the location of the position of the QFS peak, which is purely 
determined by kinematics. Therefore, we should expect that the relativistic kinematics will infiuence our results, 
especially at energies larger than 0.5 GeV. Of course, there are other dynamical relativistic effects. For energies below 
~0.25 GeV those relativistic effects were studied in neutron-deuteron elastic scattering in Ref. (S^J. There it was 
found that the combination of relativistic effects consistently incorporated is negligible below 0.1 GeV, and manifests 
itself at 0.25 GeV mostly at large scattering angles. What happens at the energies we considered when relativistic 
effects (kinematical and dynamical ones) are incorporated is uncharted territory so far, and we want to refrain from 
speculation. Work in this direction is however underway. 
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APPENDIX A: THE ip" INTEGRATION 

According to Eqs. H2.19|l . ()2.18|l . and H2.2Q|) the (/^''-integration for fixed p, q, Xp, Xq, a;^° , q" , and x" can be written 

I{Vgo,^p)= f F[cos(^" - (^,J]G[C0S(^" - ^p)] V'- (Al) 
Jo 

where the F and G arc known functions from is and T. The substitution ip' = ip" — ipq^ leads to 



I{Vq«^Vp) = / -F[C0S^']G'[C0S(¥J' - {ipp - V3go))] V 

= -^p). (A2) 



Moreover, spHtting this integral as 

Ii.Vqa -^Pp) ^ / F[cOSip']G[cOs{ip' - {(pp - iPqo))]d(p' 



2-n 

F[cos^']G'[cos(¥3' - {tfp - ^pq,,))]dLp' (A3) 
and substituting ip' ~ 2t^ ~ ip" in the second integral, one obtains 

I{^P - ^90) - F[cosp^"] |^G[cos(^" - ((^p - p^q,))] + G[cos((p" + {p^p - ^q,))]^dp" = I{\p^p - (A4) 

Consequently, the result for the (^-integration in Eq. (|A2|) or Eq. I)A4|I does not depend on the choice of the sign in 
sin((^p — fqa) = ±-\/l — cos^ [tfp — (pqo). Only cos((y5p — (pqo) is fixed by Eq. (|2.21() . and has to be known. 

Since the integral / in Eq. (|A4() depends only on the difference of the angles {ipp — ipqg), one can choose (pqg 
arbitrarily, and thus sin pq^^ and cos ipq^, required in Eq. (|2.18|l . Moreover, the trivial identities 

cos ipp = COS ipqg COs{ipp " (Pqo) - Sm ipqo sin{ipp - iPqo), 

simpp = simpq^ cos{ipp ~ y:>q„) + cos ipq„ sin{ipp - ipq„). (A5) 

are the input for cos((/3p — p?") needed in Eq. H2.20|l . The arbitrary choice of (pqg is a good check for the numerical 
correctness of the choice of variables, and we carried out those tests. 
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Eiab[GeV] 


P 


Xp 


^10 


Xq 


II 

<1,<1 


// 

X 




(Topt [mb] 


(Tei [mb] 


atr [mb] 


O-el + O-br [mb] 


0.01 


49 


4 


4 


4 


49 


4 


4 


1913.48 


1799.08 


67.81 


1866.89 




49 


8 


8 


8 


49 


8 


8 


1886.84 


1807.50 


70.14 


1877.64 




49 


12 


12 


12 


49 


12 


12 


1904.99 


1820.77 


73.75 


1894.52 




49 


16 


16 


16 


49 


16 


16 


1903.22 


1820.46 


73.20 


1893.66 


0.1 


49 


12 


12 


12 


49 


12 


12 


335.57 


259.95 


83.10 


343.05 




49 


16 


16 


16 


49 


16 


16 


343.17 


265.83 


75.84 


341.67 




49 


23 


23 


16 


49 


16 


20 


344.34 


270.05 


76.23 


346.28 




49 


23 


23 


24 


49 


24 


20 


346.16 


272.04 


76.55 


348.59 


0.5 


49 


12 


12 


12 


49 


12 


12 


40.17 


12.05 


66.32 


78.37 




49 


16 


16 


16 


49 


16 


16 


65.62 


47.76 


32.73 


80.49 




49 


20 


20 


16 


49 


16 


20 


65.93 


47.61 


38.16 


86.22 




49 


20 


16 


20 


49 


20 


16 


85.19 


61.16 


28.84 


90.00 




49 


20 


20 


20 


49 


20 


20 


85.71 


61.30 


29.86 


91.19 




49 


24 


20 


20 


49 


20 


20 


85.72 


61.24 


30.56 


91.80 




49 


20 


20 


24 


49 


24 


20 


102.17 


64.96 


33.74 


98.70 




49 


23 


23 


24 


49 


24 


20 


110.35 


64.28 


36.42 


100.70 



TABLE I: The total elastic and break-up cross sections together with the total cross section extracted via the optical theorem 
calculated from a Malfliet-Tjon type potential at two selected energies (0.01 and 0.5 GeV) as function of the grid points. The 
double prime quantities are the integration variables. The calculations are carried out in the coordinate system in which go is 
aligned parallel to the z-axis. 



Eiat [GcV] 


alf [mb] a^^, [mb] 


^!"r [mb] al [mb] al [mb] 


c^opt [mb] a«pt [mb] 


0.003 


2561.74 2561.14 


0.0 0.0 0.0 


2562.65 2562.65 


0.01 


1820.46 1820.51 


73.20 73.55 73.13 


1903.22 1902.56 


0.1 


272.04 272.20 


76.55 75.18 75.08 


346.16 346.16 


0.5 


64.28 64.61 


36.42 36.39 35.55 


110.35 110.35 


1.0 


21.90 21.90 


23.44 23.46 23.40 


49.59 49.59 



TABLE IL The total elastic cross section, total breakup cross section and total cross section extracted via the optical theorem 
calculated in different coordinate systems at selected energies. The choice of coordinate system, i.e. which vector is aligned 
parallel to the z-axis, is indicated by the superscripts qa, q and p. 
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Eiab [GeV] 


6cm [deg] 


[mb] 


^\t' [mb] 






0.01 


0.0 


537.536 


537.536 




21.8 


420.036 


420.036 




62.1 


70.726 


70.725 




93.4 


38.289 


38.289 




151.5 


227.899 


227.899 


0.2 


0.0 


676.821 


676.821 




21.8 


148.880 


148.880 




62.1 


0.363 


0.363 




93.4 


0.223 


0.223 




151.5 


0.010 


0.010 


0.5 


0.0 


519.389 


519.389 




21.8 


16.209 


16.209 




26.3 


4.430 


4.430 




62.1 


0.088 


0.088 




93.4 


0.005 


0.005 




151.5 


0.004 


0.004 


1.0 


0.0 


3.903 X 10+^ 


3.903 X 10+^ 




21.8 


5.325 X 10"^ 


5.325 X 10"^ 




62.1 


4.072 X 10"^ 


4.072 X 10"" 




93.4 


2.678 X 10"^ 


2.678 X 10~^ 




151.5 


3.705 X lO"** 


3.703 X 10""* 



TABLE III: The elastic diflerential cross sections at different energies for selected scattering angles. The cross section labeled 
by T results from the converged solution of the integral equation Eq. (12.1911 . The column labeled T' is calculated by reinserting 
the original solution into the Faddeev equation with T' = tP + tGoPT. The calculations are based on a Malfliet-Tjon type 
potential as described in the text. 



fqo [rad] 


all [mb] ol, [mb] 


<^Ht [mb] gIj,^ [mb] 


0.0 


2561.736 2561.138 


2562.649 2562.649 


TV 

5! 


2560.885 2560.729 


2562.496 2562.496 


TV 


2561.550 2561.206 


2562.091 2562.091 



TABLE IV: The total elastic cross sections at Eiab = 3.0 MeV calculated for different values of the angle tpq^ with the + sign 
of sin(<^p — v5i3o)- The calculations are carried out in two different coordinate systems, characterized by the superscripts go and 
g, which indicate, which vector is chosen to be parallel to the z-axis. 



sign[sin((^p - (^,„)] 


a™ [mb] al, [mb] 


(^Tvt [mb] cr'pi [mb] 


+ 


2561.736 2561.138 


2562.649 2562.649 




2559.674 2559.536 


2560.091 2560.091 



TABLE V: The total elastic cross sections at Eiab ~ 3.0 MeV calculated for different sign of sm{ipp — ipqo) where (pqg = 0. The 
meaning of the superscripts is the same as in Table HVI 



Va [MeV fm] 


A* A [fm ^] 


yR[MeV fm] 


^J■R [fm ^] 


EdlMeV] 


-626.8932 


1.550 


1438.7228 


3.11 


-2.2307 



TABLE VI: The parameters and deuteron binding energy for the Malfliet-Tjon type potential of our calculation. As conversion 
factor We use units such that ftc=197.3286 MeV fm = 1. 
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^0 

i 




FIG. 1: The geometry of three vectors qo, q, and p relevant in the three-body scattering problem. The independent angle 
variables Xq, Xp, and x'^^ as defined in Section II are indicated. The dashed arrows represent the normal vectors (qo x q) and 
(P X qo). 




FIG. 2: The region of singularities of the free three-particle propagator as function of the momenta q and q" . The shaded area 
in the q — q" plane indicates the region where |a;o| < 1, i.e. the region where a pole in the a;"-integration occurs. This region is 
is enclosed by the bounding curves g+ and g_ , which contain the logarithmic singularity as function of q" as given in Eqs. 14.511 
and K^ . 
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FIG. 3: The percent error in the optical theorem as function of the average distance A, of the integration grid points g" in the 
interval (0,qmaia;) at the selected laboratory projectile energies Eiab= 0.01 GeV (top panel), Eiab= 0.1 GeV (middle panel), 
and Eiab~ 0.5 GeV (bottom panel). 
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FIG. 4; The total elastic cross section a^i (dashed line), the total break-up cross section atr (dash-dotted line) and the 
total cross section evaluated by the optical theorem aopt (solid line) given as function of the projectile laboratory energy. At 
selected energies where the calculations have been carried out the sum of the calculated total elastic and break-up cross section, 
atot = Oei + Ohr, is indicated by the open diamond. The open diamonds coincide with the solid line according to the optical 
theorem, Eg. 13.231 and the numerical values are given in Table HTl 
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FIG. 5; The elastic differential cross section at 0.2 GeV, 0.5 GeV, 0.8 GeV, and 1.0 GeV projectile energy as function of the 
laboratory scattering angle. In all cases the solid line represent the full solution of the Faddeev equation. The other lines 
represent the successive sum of different orders in the multiple scattering series, short-dashed the first order, dash-dot adds up 
to the 2nd order, long-dashed to the 3rd order, and dash-dot-dot to the 4th order contribution. 
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FIG. 6: The semi-exclusive cross section at 0.2 GeV laboratory incident energy and at 24° emission angle of the emitted 
particle. The upper panel displays the entire energy range of the emitted particle, whereas the two lower panels show only the 
low and high energies in a linear scale. The full solution of the Faddeev equation is given by the solid line in all panels. The 
contribution of the lowest orders of the multiple scattering series added up successively is given by the other curves as indicated 
in the legend of the lower left panel. 
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FIG. 7: Same as Fig. El but for an angle of 39° of the emitted particle. 
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FIG. 8: The semi-exclusive cross section at 0.5 GeV laboratory incident energy and at 24° angle of the emitted particle. The 
upper panel displays the entire energy range of the emitted particle, whereas the two lower panels show only the low and high 
energies in a linear scale. The full solution of the Faddeev equation is given by the solid line in all panels. The contribution of 
the lowest orders of the multiple scattering series added up successively is given by the other curves as indicated in the legend 
of the lower left panel. 




FIG. 9: Same as Fig.lHlbut for an angle of 36° of the emitted particle. 
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FIG. 10: The semi-exclusive cross section at 1 GeV laboratory incident energy and at 18° angle of the emitted particle. The 
upper panel displays the entire energy range of the emitted particle, whereas the two lower panels show only the low and high 
energies in a linear scale. The full solution of the Faddeev equation is given by the solid line in all panels. The contribution of 
the lowest orders of the multiple scattering series added up successively is given by the other curves as indicated in the legend 
of the lower left panel. 




FIG. 11: Same as Fig. llOl but for an angle of 30° of the emitted particle. 



29 



> 

CD 



b 



> 



T-l 

5- 

b 

CM 



0.4 



0.3 - 



0.2 - 



0.1 - 



0.0 



T 1 1 1 1 f 

MT3(lst) 

AV18(lst) y/ \ 

Bonn-B(lst) / /"\ \ 



1 



1 



1 



1 



1 



1 



0.4 



390 400 410 420 430 440 450 460 

El [MeV] 

— I 1 1 



0.3 - 



0.2 - 



0.1 - 



0.0 



-| 1 1 

1st order 

2nd order 

3rd order 

4th order 

Faddeev 




390 400 410 420 430 440 450 460 



El [MeV] 



FIG. 12: The semi-exclusive cross section at 495 MeV laboratory incident energy and at 18° angle of the emitted particle. 
The upper panel displays the first order results obtained from the realistic potentials AV18 (long dashed line) and Bonn-B 
(dashed-dotted fine) together with our calculation based on the scalar MT potential of Eq. 15.21 (dashed line). The lower 
panel displays again our first order calculation from the upper panel (dashed line), together with a successive addition of the 
next three rescattering terms. The exact solution of the Faddeev equation is given by the solid line. 
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FIG. 13; The cross section (cm.) for the semi-exclusive break-up reaction in which two particles emerge in forward direction 
with a relative energy between and 3 MeV, and one particle is detected at backward angle as function of the projectile 
laboratory energy. The result of the full Faddeev calculation is given by the solid line and is compared with calculations based 
on the lowest orders of the multiple scattering series in t added up successively as indicated in the legend. 
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FIG. 14: The elastic cross section (cm.) at backward angle as function of the projectile laboratory energy. The result of the 
full Faddeev calculation is given by the solid line and compared with calculations based on the lowest orders of the multiple 
scattering series in t added up successively as indicated in the legend. 
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